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ABSTRACT 


Neglecting  magnetic  effects  the  precession  of  the  mantle  of  the  Earth 

is  responsible  for  hydrodynamic  motions  in  the  liquid  core  because,  ( i)  the 

velocity  of  the  inner  boundary  of  the  mantle  along  its  normal  is  non -zero  unless 

the  core  boundary  is  a  perfect  sphere,  (ii)  there  is  a  viscous  coupling  between 

core  and  mantle,  (ill)  there  is  a  magnetic  coupling  between  core  and  mantle 

because  lines  of  geomagnetic  force  thread  each.  In  this  preliminary  discussion 

of  the  problem,  we  suppose  that  the  core  boundary  is  an  oblate  spheroid  of 

eccentricity  e;  it  is  shown  that  the  relative  importance  of  (i)  and  (ii) 

2 

depends  in  an  essential  way  on  the  value  of  Rj  =  we  /fl  where  w  is 
the  angular  velocity  of  rotation  and  fi(«  w)  is  the  angular  velocity  of 
precession.  Since  is  almost  certainly  large  for  the  Earth,  the  present 
discussion  is  limited  to  this  case. 

The  motion  of  the  fluid  is  found  by  considering  an  initial  value 
problem  in  which  the  axis  of  rotation  of  the  spheroid  is  impulsively  moved 
at  time  t  =  0  ;  before  that  time  this  axis  is  supposed  to  be  fixed  in  space, 
the  fluid  and  envelope  turning  about  it  as  a  solid  body.  The  solution  is 
divided  into  a  steady  motion  and  transients,  and,  by  evaluating  the  effects 
of  the  viscous  boundary  layer,  the  transients  are  shown  to  decay  with  time. 

The  steady  motion  which  remains  is  a  circulation  with  constant  vorticity  in 
planes  perpendicular  to  w,  and  £  ,  the  streamlines  being  similar  and 
similarly  situated  ellipses. 


ON  THE  MOTION  OF  A  LIQUID  IN  A  SPHEROIDAL 
CAVITY  OF  A  PRECESSING  RIGID  BODY 

* 

P.  H.  Roberts 
and 

** 

K.  Stewart  son 

I.  INTRODUCTION 

Some  years  ago,  Bondi  and  Lyttleton  (1)  examined  the  motion  of  a 
liquid  in  a  spherical  cavity  of  a  rigid  body  which  was  rotating  with  an  angular 
velocity  £  about  an  axis  through  the  centre  of  the  sphere  when  this  axis,  in 
turn,  was  processing  with  angular  velocity  £(Q  «  »)  about  an  axis  fixed  in 
space.  The  aim  of  their  investigation  was  to  throw  some  light  on  the  effect  of 
the  precession  of  the  Earth  on  the  motion  of  the  liquid  core.  The  analytical 
approach  used  by  Bondi  and  Lyttleton  was  to  suppose  that  the  fluid  is 
Newtonian  and  practically  inviscid.  Then,  in  the  first  approximation, 
viscosity  was  neglected,  and  it  was  hoped  to  correct  for  viscosity  by  a  thin 
boundary  layer  near  the  surface  of  the  cavity.  Unfortunately,  on  attempting  to 
solve  the  steady  inviscid  equations,  they  arrived  at  a  definite  contradiction 
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which  strongly  suggested  to  them  that  no  ataady  atata  motion  of  a  permanent 
character  la  poaalble  for  the  fluid. 

The  apherical  cavity  praaenta  a  special  difficulty  in  the  theory  of 
rotating  fluids  because  without  viscosity  no  motion  can  be  communicated  to 
the  fluid  by  a  rotation  of  the  boundary.  The  angular  velocity  acquired  by  the 
fluid,  even  if  the  sphere  la  not  preceasing,  la  due  entirely  to  the  generation  of 
vortlcity  in  the  boundary  layer  at  the  surface  of  the  cavity  and  one  might  expect 
it  to  take  a  time  0(a2/v)  ,  where  a  ia  the  radius  of  the  cavity  and  v  the 
kinematic  viscosity,  for  the  motion  of  the  fluid  and  sphere  to  rotate  together  as 
a  rigid  body.  Let  us  suppose  this  time  has  passed,  that  the  processing  motion 
is  set  up  at  t  a  o ,  and  that  subsequently  the  fluid  may  be  regarded  as  invlscid. 
Then  while  the  solid  body  processes  the  fluid  rotates  about  a  fixed  axis,  being 
unaffected  by  the  motion  of  its  envelope.  The  motion  of  the  fluid  is  actually 
steady  but  relative  to  the  envelope  it  appears  to  be  unsteady,  its  axis  of 
rotation  processing  with  angular  velocity  -jg  and,  when  Qt  is  small, 
this  means  that  the  velocity  perturbations  are  proportional  to  Qt .  Thus  the 
contradiction  obtained  by  Bondi  and  Lyttleton  does  not  necessarily  mean  that 
the  motion  of  the  fluid  is  not  ultimately  steady.  For  a  real  fluid,  one  might 
expect  that  the  motion  described  above  would  occur  if  Qa2/v  »  1  ,  and  if, 
on  the  other  hand,  Qa2/v  «  1  ,  the  motion  of  the  fluid  and  envelope  would 
be  virtually  a  rigid  body  precession. 

This  discussion  strongly  supports  their  view  that  the  determination  of  the 
motion  is  likely  to  prove  a  problem  of  great  difficulty,  and  consequently  we 
thought  it  desirable  to  begin  our  attack  by  studying  the  related  problem  of  the 
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oblate  spharoldal  envelope  which  includes  the  sphere  as  a  special  case  but 
proves  easier  because  the  motion  of  the  boundary  of  the  cavity  must  communicate 
itself  in  part  to  the  fluid.  This  problem  is  not  without  interest  in  the  geophysical 
context  with  which  Bondi  and  Lyttleton  were  concerned  because,  although  one 
can  say  that  the  core  of  the  Earth  is  nearly  a  sphere  in  that 


9 


where  a  and  b  are  the  semi-major  and  minor  axes  of  the  core,  one  cannot  be 
sure  of  the  values  of 


V 


Oa 


bfh, 

2 
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which,  with  R3  =  wa  /v  (  »  1)  ,  control  the  motion  of  the  fluid  in  the  core 
relative  to  the  rigid  body  rotation.  If  R^  =  0,  the  motion  is  critically 
dependent  on  R2  while  if  R^  is  large  the  motion  is  independent  of 

except,  possibly,  if  rir2  ~  1  •  Our  investigation  here  is  based  on  the 

assumptions  that  R ^  »  1 ,  Rj  R  2  »  1  ;  it  is  hoped  to  explore  other  limiting 

situations  in  a  later  paper.  These  assumptions  are  not  unreasonable  in  the 

2  2  2 

geophysical  context  since,  for  the  surface  of  the  Earth,  (a  -b  )/a  «  1/300  , 

and  the  value  of  this  parameter  for  the  core  is  not  likely  to  be  an  order  of 

magnitude  less.  Further,  w/ft  ~  5. 102 * * * 6 * 8  ,  and  v  is  almost  certainly  much 

8  2 

less  than  10  cm  /sec.  However,  it  should  not  be  forgotten  that  our  results 


2  2  Z 

are  extremely  sensitive  to  (a  -  b  )  /a  when  a  »  b  .  It  may  also  be 
appropriate  to  remark  that,  as  a  geophysical  model,  our  present  approach 
Involves  a  gross  oversimplification  of  the  dynamics  of  the  core  since  it 
excludes  magnetic  forces,  and  these  are  likely  to  be  the  most  potent  agency 
in  coupling  core  and  mantle  (cf. ,  e.  g. ,  Bullard,  2).  Of  course,  a  realistic 
model  including  magnetic  forces  would  be  immensely  complicated  (although, 
in  its  favour,  it  should  be  mentioned  that  the  case  a  «*  b  would  no  longer 
be  singular!).  Nevertheless,  we  believe  our  simple  model  has  some  bearing 
on  the  generation  of  the  geomagnetic  field.  It  has  long  been  known  ( 2)  that 
the  only  processes  which  are  energetically  capable  of  generating  the  geomagnetic 
field  by  dynamo  action  are  thermal  convection  and  precession.  Bullard  ( 2) , 
while  himself  favouring  the  former  of  these  mechanisms,  was  unable  to  rule 
out  the  latter.  Indeed,  the  subsequent  suggestion  (1)  that  turbulence  might 
be  generated  by  the  processional  effects  near  the  singularities  of  the  boundary 
layer  at  latitudes  30*  N  and  30*  S  made  it  seem  (3)  not  unlikely  that  pro¬ 
cessional  motion  played  an  Important  part  in  the  geomagnetic  dynamo.  However, 
turbulence  was  invoked  by  Bondi  and  Lyttleton  as  a  motion  which  the  fluid  would 
be  likely  to  adopt  in  the  absence  of  a  permanent  steady >  state  flow.  In  our  model, 
on  the  other  hand,  a  steady  state  flow  does  exist,  and  turbulence  could  arise  only 
if  it  is  unstable  ( see  below) .  If  we  exclude  this  possibility,  it  is  clear  that, 
even  though  in  our  model  the  order  of  magnitude  of  the  pieces  sional  motions 
( ~  1  cm.  / sec.  near  the  core  boundary)  should  suffice  to  make  their  inductive 
effects  felt,  their  structure  is  far  too  simple,  in  regard  to  the  restrictions  imposed 
by  Cowling's  theorem  (4) ,  to  maintain  dynamo  action.  Other  motions  would  be 
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requlred  to  complete  the  regenerative  oyoU.  In  principle,  even  a  simple 
differential  rotation  about  tha  axis  of  £  might  suffice  (5),  and  It  may  ba 
worth  noting  that  such  motions  might  arlsa  If  «  1  in  tha  core. 

Evan  for  an  oblate  spheroidal  cavity  containing  almost  lnvlscid  fluid 
tha  determination  of  the  flow  has  some  unusual  features.  Thus  the  governing 
equations  for  an  lnvlscid  fluid  are  hyperbolic  while  tha  boundary  condition  is 
a  relation  between  the  function  and  its  normal  derivative.  Such  a  boundary 
condition  is  appropriate  to  an  elliptic  differential  equation  and  there  are  no 
corresponding  existence  and  uniqueness  theorems  for  hyperbolic  equations. 

In  fact  one  can  construct  examples  of  non-existence  and  non -uniqueness  as 
follows. 

Suppose  that  the  governing  equation  is 

8x  0y 

2  2 

and  that  (a)  #*0  on  the  circle  x  +  y  si,  or  (b)  on  three  sides 

of  the  square  |x|  <  1,  lyl  <  1  and  #  =  1  on  the  fourth  side.  One  solution 
of  (a)  is 

+  =  A(x2  +  y2  -  1) 

Inside  the  circle,  where  A  is  arbitrary.  This  solution  can  be  added  to  any 
other,  so  that  there  is  no  unique  solution  of  (1.  1)  in  this  case.  Again  every 
solution  of  (1.  1)  which  vanishes  on  three  sides  of  a  square  must  vanish  on 
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the  fourth  aid*  ao  that  tha  aolution  of  problem  (b)  cannot  axiat. 

Fortunately  In  the  preaent  problem  there  la  a  simple  physical  explanation 
of  these  phenomena.  The  non-uniqueness  is  associated  with  the  existence  of 
free  oscillatory  motions  of  the  fluid  in  which  the  boundary  remains  fixed  and 
the  non-existence  is  associated  with  the  resonance  which  can  occur  when  the 
boundary  is  made  to  oscillate  with  the  same  period  as  that  of  one  of  the  free 
oscillations  of  the  fluid. 

The  difficulty  about  the  inappropriate  type  of  boundary  condition  is 
removed  by  considering  motions  which  start  from  a  position  of  relative  rest. 

We  suppose  that  the  fluid  and  envelope  are  initially  rotating  as  if  solid  with 
angular  velocity  £>  and  that  at  time  t  =  0  the  envelope  is  set  processing  with 
angular  velocity  Jg .  The  equations  governing  the  induced  motion  of  the  fluid 
are  unsteady  but  the  time  dependence  is  removed  on  applying  a  Laplace  trans¬ 
formation  with  parameter  s  and  the  equation  is  reduced  to  a  single  one  of 
Laplace's  form  apart  from  a  scaling  factor  which  is  a  function  of  s.  A  general 
solution  can  be  written  down  and,  on  inversion,  the  history  of  the  motion  can  be 
traced.  From  an  examination  of  the  solution  at  large  times  we  can  relate  the  non- 
uniqueness  to  free  oscillations  set  up  by  the  initial  motion  and  the  non-existence 
to  resonance. 

For  a  spheroid  the  actual  inviscid  solution  is  straightforward  being 
effectively  given  by  Hough  and  by  Poincarf  (6),  the  difficulty  found  by  B<~  ol 
and  Lyttleton  for  a  *  b  is  seen  to  be  a  resonance  and  explicable  on  the  ...  . 
stated  in  physical  terms  at  the  beginning  of  this  section. 
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ln  |4  the  boundary  layar  dua  to  this  motion  is  investigated.  Tha 
lnvlsold  motion  outside  consists  of  two  parts  (i)  dspandant  on  tha  aslmuth 
angle  6  and  dua  to  tha  steady  motion  of  the  boundary,  (ii)  independent  of 
8  and  dua  to  tha  angular  velocity  which  the  boundary  acquired  Initially  but 
failed  to  communicate  Immediately  to  the  fluid.  The  boundary  layer  associated 
with  ( i)  was  studied  by  Bondi  and  Lyttleton  but  we  repeat  their  investigation 
here  since  we  are  including  time  as  a  parameter.  It  develops  a  singularity  on 
certain  circles  of  latitude  as  s  approaches  any  value  on  the  imaginary  axis 
between  >31  and  i.  The  boundary  layer  associated  with  (ii)  has  also 
been  studied  earlier  by  Stewartson  (7)  and  Proudraan  (8) ,  and  it  develops 
singularities  as  s  approaches  any  value  on  the  imaginary  axis  between  ±21 . 

Bondi  and  Lyttleton  (1)  noticed  that  the  normal  velocity  on  the  critical 
circle  becomes  infinite  in  the  limit  and  suggested  that  the  flow  near  it  was 
unstable  leading  to  turbulence  here  and  elsewhere.  Although  one  cannot  be 
sure  without  carrying  out  an  actual  experiment,  more  recent  investigations 
do  not  support  this  view.  The  singularity  on  a  critical  line  is  a  singularity 
in  the  boundary  layer  sense,  i.  e. ,  if  vn  is  the  normal  velocity  just  outside 
the  boundary  layer,  the  equations  imply  that  vr  ■  O(v^)  almost  everywhere 
and  on  approaching  the  critical  line  v"^  vn  -*  « .  This  does  not  necessarily 
mean  that  vr  -*  »  and  indeed,  in  a  related  problem  (9),  a  more  detailed 

4- 

investigation  of  the  flow  properties  near  the  critical  line  implied  v  =  O(v  ) 

n 

which  is  still  small  in  the  limit  v  -*  0 .  Further  a  detailed  solution  of  a 
related  problem  to  (11)  by  Stewartson  (10)  elucidates  the  role  of  the 
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critical  circles  as  the  origin  of  shear  layers  which  penetrate  into  the  fluid, 
transporting  fluid  from  one  part  of  the  boundary  to  another  via  the  interior 
and  also  adjusting  the  angular  velocity  of  the  main  body  of  fluid  when 
required.  Free  shear  layers  are  certainly  notoriously  prone  to  instability; 
however,  after  Proudman  (8)  had  postulated  their  existence  in  e  discussion 
of  the  flow  between  two  concentric  spheres,  experiments  carried  out  to  tost 
his  conclusions  did  not,  in  all  cases,  reveal  Instability.  It  was  found  (11) 
that,  if  the  angular  velocity  of  the  inner  sphere  exceeded  that  of  the  outer 
sphere,  the  shear  layer  was  rapidly  destroyed  by  Taylor- Gortle  instabilities 
and  by  the  formation  of  Kdrmdn  vortex  sheets  at  the  layer.  If  however,  the 
outer  sphere  rotated  more  rapidly  than  the  inner  sphere,  the  shear  layer  showed 
scarcely  any  tendency  to  instability. 

In  §5  and  §6  the  (tertiary)  modification  to  the  secondary  flow 
induced  by  these  boundary  layers  is  discussed  with  particular  reference  to 
the  oscillations  produced  by  the  initial  motion  of  the  envelope.  It  is  shown 
that,  of  the  motions  dependent  on  8 ,  only  the  steady  motion  survives  as 
t  -*  00 .  The  effect  of  the  boundary  layer  Independent  of  8  is  to  lead  to  a 
breakdown  in  the  tertiary  flow  as  t  and  it  is  argued  that  this  must 
mean  that,  relative  to  the  boundary,  the  motion  independent  of  8  must  die 
out  as  t  -*  «  due  to  the  communication  of  vorticity  to  the  fluid  via  the 
boundary  layer. 
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2.  EQUATIONS  OF  UNSTEADY  MOTION 

We  consider  a  mass  of  incompressible  fluid  which  occupies  the  whole 
of  a  rigid  envelope  whose  internal  boundary,  SQ  ,  is  a  surface  of  revolution 
with  axis  of  symmetry  Lg  .  Initially  the  fluid  and  the  envelope  are  rotating 
about  the  axis  with  angular  velocity  £  which  is  an  absolute  constant.  At 
time  t  =  0 ,  the  axis  Lg  is  set  rotating  with  a  small  uniform  angular  velocity 
£  about  an  axis  Lg  fixed  in  space  which  Intersects  Lg  at  a  point  O,  and 
which  is  inclined  at  an  angle  a  to  it.  For  definiteness,  we  shall  suppose  that 
the  perturbed  motion  is  started  impulsively  so  that  £  is  also  an  absolute 
constant.  We  wish  to  find  the  subsequent  motion  of  the  fluid. 

Consider  a  reference  frame,  ,  rotating  with  angular  velocity  £ 
relative  to  which  Lg  and  Lg  are  stationary,  and  let  the  velocity  of  the 
fluid  relative  to  be 


£  =  •  (2.1) 

At  t  =  0- ,  J4,  »  0 ,  and  at  t  =  0+ ,  #  0  in  virtue  of  the  impulsive  motion 

of  Sq  ,  but  it  will  be  determinate  from  the  initial  motion  of  SQ  .  The  actual 
velocity,  x  >  of  the  fluid,  relative  to  axes  fixed  in  space  and  Instantaneously 
coinciding  with  ,  is 


(2.  2) 
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We  know  that,  at  t  *  0+ ,  ^  +  Jg  1*  irrotational  and  uniquely  determinate 

from  the  motion  of  SQ ,  and  hanoa  la  alto  known  at  t  ■  0+  .  Tha  aquations 
of  motion  ralatlva  to  ara 

x  X  +  (&•  Of«d)X  *  -orad(  W  +  £)  +  vV2£  ,  (2.3) 

dlv  x  a  0  *  ( 2. 4) 

where  p  Is  the  pressure,  p  the  density,  v  the  kinematic  viscosity,  and 
W  the  potential  per  unit  mass  of  the  external  forces,  supposed  conservative, 
which  act  on  the  fluid.  We  now  express  (  2.  3)  in  terms  of  ,  obtaining 


+  *(£  +  &)*  Jt  *  U  +  £  X,*)  X  curl  A 


=  -grad 


w+£  -  y {<>!»  +  je>x^t>2 
L  +  (&xx,)  + 


+  vV  Ji  +  ,£  X(£  Xjsj»  ) 


whence,  on  neglecting  squares  and  products  of  j£,  end  jg  ,  we  find 


+  2£Xji  -  (£X,£)X  curl  &  »-«  grad  V'  ♦  vV2%  -  ,  (2.5) 

where 


«V«  »  W  +  £  +  ±£2  -  (£X,£)-2  +  ,l£)r2  • 


(2.6) 


#374 


-11- 


This  equation  dllfara  from  that  given  by  Bondi  and  Lyttleton  (1),  viz. 

+  2^X  X^)  Xeurt^  = -w  grad  v  +  ,  (2.7) 

where 

«V  -  W  +  ^+  -  (*>X,£)2  +  2(£-£)(Q'£)  &)*  • 

The  first  term  In  ( 2. 7)  In  brackets  was  omitted  by  Bondi  and  Lyttleton 
because  they  assumed  that  the  motion  is  steady.  So  far  as  the  motion  of  the 
fluid  is  concerned,  the  only  significant  difference  between  ( 2.  5)  and  ( 2. 7) 

Is  the  factor  2  in  the  last  term  of  ( 2. 5)  .  We  shall  show  ( $3)  that  this 
factor  is  crucial  for  obtaining  a  consistent  solution  of  the  problem  when  SQ 
is  a  sphere.  The  equation  of  continuity  ( 2. 4)  is  equivalent  to 

divji"0  •  (2.8) 

Now  choose  a  system  of  cylindrical  polar  coordinates  (r,  6,  z)  in 
&  with  Lg  as  the  z-axis,  r  measuring  the  distance  of  a  typical  point 
P(^£)  from  Lg,  and  8  being  the  angle  between  the  plane  defined  by  P  and 
Lg  and  the  plane  defined  by  Lg  and  Lg  .  Then 

18  •/■£)£  =  (°»  °»  COS  a  -  n«r  sin  a  cos  8)  ,  (2.9) 

where  O  =  -  |jQ,|  is  the  (retrograde)  angular  velocity  of  precession. 
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To  begin  with,  neglect  viscosity,  and  danota  tha  velocity  by 
with  component*  (uof,  uQe,  uQ|5 )  raapaotivaly  along  tha  directions  of 
r,  6  and  x  increasing;  (2.5)  and  (2.8)  are  than  equivalent  to 


8u__  8u 


8t 

-f  co 

T  88 

8un4 

08 

8u„ 

8t 

T  —• 

88 

8uft 

Oz 

.  ,  %. 

8V" 

1  S  —  (ilam- 

*08  dr 


«  8V» 

*0r  “  “F  88  * 


8V" 

“  “  “"IF  +  2w0r  8in  “  co*  ®  » 


(  2. 10) 


1  8  . 

Feft'V 


1  %8 

88 


8u 


Oz 

8s 


where 


V"  *  V*  -  Ox2  cos  a  -  run.  . 

Wv 

The  boundary  conditions  require  that  the  normal  component  of  the  fluid  velocity 
is  xero  on  SQ ,  and  also  that,  at  t  =  0+ ,  )t  +  X,  is  lrrotational  with 
velocity  potential  (say). 
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In  order  to  solve  ( 2. 10)  It  Is  convenient  to  take  the  Laplace  transform 
of  the  dependent  variables  with  respect  to  time, denoting  the  result  by  a  star, 
e.g., 

00 

u*r  ( r,  e,  z,  s)  ■  /  e““*tu0r(r,  ©,  z,  t)dt  .  (2*11) 


The  governing  equations  then  reduce  to 


1  8V»* 
r  80 


+  w  tlWi  ’ 


(2.12) 


:(U0z>i  + 


2 Osin  a 

cos 


rcos  0 


t 


where  the  suffix  1  denotes  initial  conditions.  Since,  at  t  =  0+ , 

+  &x,i  is  derivable  from  the  potential  ^  ,  equations  ( 2. 12) 
may  be  simplified,  on  writing 

V*  =  VM*  -  —  +  —  rz  sin  a  sin  8  (2. 13) 

w  Ti  w 


to 
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where  uQr,  uoe,  u0c,  and  V  ara  independent  of  0.  On  substituting  into 
(2.14)  and  (2.15),  we  obtain 


n 


#374 


Or 


M  (•  +1)  *  +  4 
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(2.17) 


-  21  Osin  a _ 1_  8V 

U0s  “  "  us  r  “  (s  +  i)  8z 


where  V  satisfies 


1__8_ 
r  8r 


.v  Ui+iii+ii  a*  =0 

r  (s  +  i)  8* 


(2. 18) 


If  SQ  Is  given  by  f(r,  z)  si,  the  boundary  condition  to  be  satisfied  is 


-  8f  —  8f 
U0r  8r  +  U0z  8z 


s  0 


at  f  si 


(2.19) 


I 
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3.  THE  SEOONDMff  INVX8CXD  FLOW  FOR  A  SPHEROID 

The  problem  posed  by  equations  (2. 17)  to  (2. 19)  can  be  reduced  to  a 
solution  of  Laplace's  equation  on  writing  [(s  +  i)  +4]/(s  +  i)  ■  %  ,  and 
replacing  z  by  £z  .  We  can  expe  ct,  therefore,  that  it  can  always  be  solved 
for  the  boundary  condition  ( 2. 19),  which  is  equivalent  to  a  linear  relation 
between  V  ,  its  normal  derivative,  and  possibly  also  its  tangential  derivative. 
Having  obtained  the  solution,  the  ultimate  flow  may  be  found  by  letting 
s  -*  0  ,  and  taking  note  of  the  poles  and  branch  points  in  the  solution  in 
the  half  plane  (R.{  s)  >  0  . 

In  this  paper  we  are  particularly  concerned  with  the  properties  of  the 
solution  when  the  envelope  SQ  is  the  oblate  spheroid 


^  \  -1,  («>b)  ,  (3.1) 

a  b 

of  which  Bondi  and  Lyttleton  considered  the  special  case  t>b,  For  the 
boundary  (3.1),  (2.19)  becomes 


"7  “or  +  f?  “o.  ”  0  ’ 

ft  D 


and  the  appropriate  solution  of  ( 2. 18)  is  found  by  writing 


(3. 2) 


V  =  Arz 


f 


(3.3) 
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According  to  ( 3.  6)  to  (3.8)  ,  the  secondary  motion  of  the  fluid 
consists  of  two  parts:  one  is  essentially  due  to  the  initial  motion  of  the 
boundaries;  the  other  is  steady.  The  first  part  can  itself  be  divided  into 
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two  components.  One  is  rigid  body  rotation  of  angular  velocity  O  cos  a 
about  Lg  which  the  rotation  of  the  boundary  failed  to  communicate  to  the 
fluid.  The  other  is  a  free  oscillation  of  the  fluid.  The  steady  part  of  the 
motion  consists  of  the  rigid  body  rotation  ju>  (cf.  eq.  [2.1]),  and  a 
circulation  ( given  by  the  terms  in  sin  6 ,  cos  8  and  sin  8  in  [  3.  6] , 

[3.7]  and  [3. 8],  respectively)  in  planes  perpendicular  to  Lg  the 
streamlines  being  similar  and  similarly  situated  ellipses.  The  vorticlty 
of  this  circulation  is  constant,  and  equal  to 

-  20  sin  a  ^  +  b  ^  . 

(a  -b  ) 

In  the  case  of  a  sphere  a  *  b ,  the  relevant  periods  of  free  oscillation 
of  the  fluid  are  infinite,  so  that  a  resonance  develops  and 


uQr  =  -  Oz  sin  a( sin  8  +  wt  cos  8) , 


uQ8  =  Oz  sin  a(cos  6  +  «t  sin  a)  +  Or  cos  a  , 


\  (3.9) 


uQz  =  Or  sin  a(sin  8  +  wt  cos  8)  ,  J 

or,  in  vector  notation, 


Ji  =  -£x,L  -  M£x  «) 


(3.10) 
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The  physical  Interpretation  of  this  rasult  la  that,  as  on*  sxpocts,  tha  motion 
of  a  spherical  boundary  exerts  no  Influence  on  the  motion  of  the  fluid  and 
hence,  relative  to  the  frame  $  ,  the  axis  of  rotation  of  the  fluid  rotates 
with  angular  velocity  •£  about  Lg  .  Consequently,  after  time  t(Qt«l), 
the  fluid  rotates  about  an  axis  making  an  angle  of 

*«*«  | 

1*1  I 

with  Lfi  ,  and  this  is  indicated  by  ( 3. 10) .  It  is  noted  that,  on  substituting 
(3.10)  Into  (2.7),  the  equations  of  motion  are  satisfied  identically,  and  so 
any  change  In  the  motion  of  the  fluid  must  be  Initiated  through  a  boundary 
layer  arising  from  the  non-satisfaction  of  the  requirement  that  ^  »  0  on  the 
boundary  SQ  .  Had  the  equation  ( 2. 7)  required  by  Bondi  and  Lyttleton  been 
used,  t  would  be  replaced  by  t  In  (3.10),  and  the  new  formula  for 
would  not  be  consistent  with  this  physical  argument. 
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4.  THE  BOUNDARY  LAYER 

The  inviscid  solution  obtained  in  ths  prsvious  •  Action  satisfies  not 
only  the  inviscid  equations  but  also  the  viscous  equations  (2. 5).  Viscosity 
therefore  manifests  itself  only  through  the  fact  that  this  solution  does  not 
satisfy  the  viscous  boundary  condition  ^>0  on  the  spheroidal  envelope 
SQ  ,  and  leads  to  a  boundary  layer  whose  thickness  we  can  anticipate  to 
be  O(v^)  if  v  is  small  where  the  adjustments  in  the  tangential  components 
of  velocity  are  made.  We  shall  calculate  this  boundary  layer  in  the  present 
section  and  its  effect  on  the  inviscid  flow  in  the  interior  will  be  discussed 
in  the  following  sections. 

Let  us  write 

A*Ao  +  £i’  (4*l> 

where  is  defined  by  (2. 1) ,  jj,0  is  its  value  if  v  »  0  ,  and 
is  the  correction  due  to  viscosity.  Then,  to  order  v  *  , 
is  the  contribution  from  the  boundary  layer  to  ^ .  Introduce  spheroidal 
polar  coordinates  (X^,  |ij,  8),  where  and  ^  are  defined  by 

r»(\J  +  c2)*(l-i*J)*  ,  (4.2) 

c2=a2-  b2  and  the  spheroid  8Q  is  given  by  *  b  .  Then,  using  (3. 6)  to 
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(3. 8),  we  see  that  the  boundary  condition  =  0  requires  that,  at 
“  b  , 


20ab  sin  a 

2  2 
(a  -b‘) 


[sinG - - — =-  sin(8  +  wkt)](b2  +  |i2c2)^ 

(a  +b  ) 


(4.3) 


208^ 


sin  a 


[cos  6-  — = - T-  cos(0  +  wkt)  ](i,  -  Oa  cosafl-iO5 

(a2  +  b2)  1  1 


Beyond  the  boundary  layer,  when  (b  -  Xj)v^  is  large  and  positive,  =  o(l)  j 

we  shall  then  be  particularly  interested  in  u,  which  is  O(v^)  and,  unlike 

1 

the  other  components  of  velocity,  is  not  exponentially  small.  Consequently 
it  engenders  a  tertiary  flow  through  the  interior  of  SQ  .  From  the  boundary 
condition  (4.  3)  the  flow  in  the  boundary  layer  may  be  divided  into  two  parts 
one  of  which  (^2)  is  dependent  on  6  ,  and  one  which  (j^)  is  independent 
of  6 .  We  consider  the  part  dependent  on  8  first. 

Write  the  Laplace  transform  of  ^  with  respect  to  «t  as 

(JI2  e10)  ,  (4.4) 

and  in  the  boundary  layer  make  the  conventional  assumption  tw<u  iu«=  opera. or 
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8/8X^  *  O(v^)  while  8/8|ij  and  8/88  art  0(1).  Than  as  it  usual  with 
Ekman  flows,  in  tha  governing  aquations  6  and  ^  may  batakan  constant. 

Xj  may  ba  replaced  by  b  except  whan  it  appears  in  8/8Xj .  Tha  governing 
aquations  reduce  to 


_  V  _  a2v  '  *>1 

(•+1)^2^  +  2  2  2  4  U28  *  2  2  2  2 

^1  (b  +c  |ip*  «(b%c|ip  8X^ 


(4.5) 


(s  +  i)u2e 


(b2tc2,^  *1 


2 

a  v 


2  2  2 
«(b%c  lip 


(4.6) 


a(b  +c  |i 


8u2x 

2.  2  2,i  *1  .  8  2^2.  2  2.i-  ,  . 

+e  <*•>  —  +  <*»  +e  1*1 '  + 


—  2  2  2 
iu2e(bSc  lip 

•0-1*?)* 


*  0  (4.7) 


and  tha  boundary  conditions  to 


U  Us 

U2Xj  =  0,  U2|i 


2iQab  sin  a(b2+c2|*?)*(s  +  i)  _  208^  sina(s  +  l)|i^ 


1  «s[s(a2  +  b2)  -  i(a2-b2)] 


20  «s(s(a2  +  b2)-i(a2-b2)] 


at  X,  =  b.  Also  u_  and  u,fl  *•  0  as  (b-X,)v 


i 


#374 


-23- 


Integrating  (4.5)  (4.6)  we  have 


»  .  o«bl  .in  ,<.^1)  J[.vV+cV?«.li2l>ll'bl 

*1  «s[s(a2+b2)-i(a2-b2)]  1  1 


ML.  b) 

-  (an2-  N/(b2+c2n2)]a  1 


u29*  rrm  it:  W'tAfc 


«s[ s(c  +b  )  -i(  a  -  b  )] 


2.  2  2.,  62(Xl_b) 


2.  2  m  *i*xrb| 


+  (atij-N/jb  +c  jij)]® 


(4.8) 


where 


2  u  2  2  2  .  2***ia 

‘2 - j(b  ♦«  !.,)[•  +1  -  ,  ;  ;  2  ] 

va4  ^(b  +  c  m>1  ) 


.2  »  ,t2  .  2  2„  .  ,  .  ^1 

61  “ — 2  (b  +c  M-i>£  8  +l+  “ — l - — 

va4  V(d  +  cn4) 


]  , 


(4.9) 


and  the  signs  of  the  square  root  of  the  expression  for  6^  and  62  are  decided 

by  requiring  that,  when  s  is  large,  real  and  positive,  6^  and  62  are  both 

positive.  The  value  of  p,  now  follows  by  the  direct  integration  of  (4.7) 

ZX1 

so  that  all  properties  of  the  boundary  layer  are  now  formally  known.  It  is  noted 
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that  tha  assumptions  on  which  this  boundary  layer  theory  is  built  are  consistent 

provided  only  that  6^  and  62  do  not  vanish.  Exceptional  cases  arise  only  if 

s  is  purely  imaginary  and  Is  |  <  3  ,  when  the  assumptions  break  down  on 

certain  circles  ^  *  constant.  Conceivably  this  may  be  serious  at  large  times 

since  singularities  occur  in  the  s -plane  on  the  imaginary  axis}  on  inverting 

we  could  get  contributions  to  of  order  tne  ,  (n>0).  However, 

it  is  believed  that  this  does  not  happen.  In  a  closely  related  problem  ( 9) 

the  flow  has  been  studied  in  the  neighbourhood  of  these  lines,  and  it  has 

been  shown  that  main  effect  is  to  change  the  thickness  of  the  layer  from 

0(  v*)  to  0(  v*) .  Adapting  the  argument  to  the  present  problem  similar 

results  are  obtained,  and  it  is  found  that,  in  the  neighbourhood  of  the 

4*  4 

singular  lines,  =  0(  v* ) ,  instead  of  0(  v*)  as  elsewhere  in  the 

"  i 

boundary  layer.  Further,  on  computing  the  contribution  from  the  neighbourhood 
of  these  lines  to  the  flow  in  the  deep  interior  of  the  spheroid,  we  find  it  to  be 
0(  v+)  as  against  0(  v^)  from  the  rest  of  the  boundary  layer. 

Of  particular  interest  is  the  value  u_  just  outside  the  boundary 
layer  1.  e.  when  (b-Xj)v*  is  large,  but  (b-X^)  is  small.  This  is  obtained 
by  integrating  (4. 7)  with  respect  to  Xj  from  -»  to  b  and  using  (4. 9). 
Denoting  the  value  of  u^  by  u^  (b,  6)  we  have 
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,  2  2 
+  C 


(b,  n.,  6)  = 
1  1 


Qabl  sin  a(s  +  i) 


»a[s(a2+b2)  -  l(a2-  b2)J 


2.  2  2. 


8 

8^1 


(b  +c  up 
a(l-n2)*  ^ 


<l-^(b2+cV 


2A  2  2*i 

%  J.  /«  II  1  • 


a^+N/^+c2**2)  ajij-^b^+c^j) 

—  r 


2.  2  2. 


2  .  2  2i\ 

A.  r*  ii  l  •  I 


a*t,  +(b  +  c  VTr  aj*.  -  (b  +c  |i.“) 

+ 


i 


(4. 10) 


Now  let  us  consider  the  second  part  ( Ji3)  of  the  boundary  layer  due  to 

that  part  of  (4.3)  independent  of  6.  Let  the  Laplace  transform  of  ^  with 

e 

respect  to  ut  be  JI3  •  Then 


u 


3X, 


0, 


u 


3*, 


1  2  i 

=  0,  uje  *  -  Qa  cos  atl-Pj)*  /  us 


(4.11) 


and  by  an  analogous  argument  to  that  given  above  the  equations  governing  the 
flow  in  the  boundary  layer  reduce  to 


*u  u  +  2  2  2i 

(bSc  Vp* 


U30  = 


2 

a  v 


3**, 


w(b^  +  c2»i2)  8X* 


(4. 12) 


su 


3l- 


Z^a 


*  - 


(b2  +  c2»*p 


7? 


2 

a  v 


2  2  2  2 
«(b‘  +  c  pp  8X‘ 


(4. 13) 
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8U3* 


•<b2+cV^ [  a  - •‘f )*  (b2+cVf)1u3»  j .  o 


(4.14) 


Integrating  (4.12)  and  (4.13)  we  have 


u35j  “  "  2us 


Qalcosa(l-nf)*  6  (X.  -b)  6.(X.-b) 

- i—  [e  3  1  -  e  4  1  ] 


nacosa(l-iij)*  6  (X  -b)  64{X  -b) 

u3$  - - SS -  <e  +  6  1 


(4. 15) 


where 


e2  «  ..2,  2  2.  . 

&3  =  — ^  (b  +  c  HjMs  - 

va 


Zi^a 

(b2  +  c2»*2)^ 


.2  «  .,2,  2  2,, 

®4  ““2  (b  +C  V*8  + 
va 


zi^a 

(b2  +  c2^ 


(4.16) 


and  6^  and  6^  are  real  and  positive  when  s  is  real,  positive  and  large. 

The  consistency  of  this  solution  is  subject  to  the  same  qualifications  as 

(4.10).  Of  particular  Interest  here  too  is  the  value  of  u,*  just  outside 

*X1 

the  boundary  layer,  and  denoting  it  by  u3  *  (b,  Hjj  8)  we  have 
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5.  THE  TERTIARY  INVISCID  FLOW 

The  boundary  layer  solution  obtained  In  §4  shows  that  there  must  be 
a  normal  velocity  at  the  edge  of  the  layer  (1.  e.  as  •*  b-  on  the  inviscid 
scale)  of  order  and  given  by  the  sum  of  (4.10)  and  (4.17).  In  turn 
this  must  induce  a  tertiary  inviscid  flow  throughout  the  interior  of  the  spheroid, 
i.  e.  a  flow  governed  by  ( 2. 14)  and  ( 2.  IS)  excluding  the  forcing  terms,  and 
satisfying  (  2. 19)  with  a  right  hand  side  effectively  equal  to  the  sum  of  ( 4. 10) 
and  (4.17)  (instead  of  being  zero).  Since  (4.10)  is  dependent  on  6  while 
(4. 17)  is  not,  it  is  convenient  to  treat  their  contributions  separately  and  in 
this  section  we  shall  consider  the  consequences  of  (4.10).  In  §6  below 
we  shall  consider  the  consequences  of  (4. 17) . 

Denoting  the  corresponding  velocity  in  the  spheroid  by  jj,  we  have, 
from  (2.17),  that  ^  can  be  expressed  in  terms  of  a  scalar  V4  which 
satisfies  (  2. 18),  the  forcing  term  in  ( 2. 17)  being  again  set  equal  to  zero. 
Hence,  after  stretching  the  z  coordinate,  V 4  satisfies  Laplace's  equation 
and  the  appropriate  solution  can  be  formally  written  down.  Let 


r  = 


[  (  s  +  i)  +  4] 


1 


(1 


z  = 


(s  +  i)  **4X4 


(5.1) 


where 


2 

*4 


2  12  2 
a*  +  ^c  (s  +  i) 


(5.  2) 


#374 


-29- 


In  terms  of  and  X  4  tha  spheroid  is  given  by 

X4«i{*  +  i)b  (5.3) 

and  on  the  spheroid  =  K  •  The  most  general  acceptable  form  for  V4 

is 


u-^)* 


(*4  +  v2)*  E  a^pk~)p:(^) 


co 


iX 


n=l 


n  n' 


n' 


(5.4) 


the  P  being  Legendre  polynomials  and  the  A  constants  to  be  found, 
n  n 

Further,  the  boundary  condition  associated  with  (4. 10),  viz. 


r 

2 

a 


v+ 


2  u4z 


2  2  , 

A  D 


(b,  *,  0) 

1 


f 


becomes,  in  terms  of  , 


a8V4 


2ib 


2  0X  '  2  '  “  V4  =  (b2  +  cV)*  u 

2  8X4  a[(s  +  i)  +  4  4  4X1 


at  X A  =  (s  +  i)b/2  .  It  follows  from  (4.10)  that  u.. 

4  4XX 


(5.5) 


is  an  odd  function 
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2  4 

of  >  and  behave  a  like  (1  -  u  )*  near  n  ■  ±  1 .  Consequently  n  is  even 
In  { 5. 4) ,  and  the  various  An  are  determined  from  the  equation 

«Bn  ,  (5.6) 

with  x4=(8  +  i)b/2  and 


V 


a  _9_ 
2  8X  , 


2  2  4  ^4 

i(k4+vp;‘-7±'i ♦ 


2ib(X  +y)a  iX 

- t-TL-*nlT±) 

a[(s  +  i)  +4]  Y4 


2n(n  + 1) 
2n+ 1 


+1  , 

= /  (b2+ 
-1 


2  2.j  - 
c  H  )*  u 


dH 


(5.7) 


From  this  point  on,  the  determination  of  is  in  principle  straight¬ 
forward  but  involves  the  evaluation  of  complicated  inversion  and  other  integrals. 
However,  since  v  is  assumed  to  be  small  the  effect  of  the  tertiary  flow  may 
be  neglected  unless  the  Ar  have  poles,  qua  functions  of  s,  in  the  half 
plane  <&  (s)  >  0  ,  or  sufficiently  strong  singularities  on  the  imaginary  axis 
of  s  .  Our  purpose  in  the  rest  of  this  section  is  to  give  arguments  to  exclude 
such  possibilities. 

First  we  show  that  A  has  no  singularities  in  the  half  plane 

n 

&  (s)  >  0  .  A  singularity  can  arise  if  either  the  coefficient  of  A  in 

n 

(5.6)  vanishes  or  if  Bn  is  singular,  or  if  (2.18)  breaks  down.  Taking 
these  possibilities  in  order,  the  coefficient  of  A^  vanishes  if 
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<s  -  l)bP^(^)  +  2iY4n(n  +  l)Pn(*)  =  0, 
using  the  properties  of  Legendre  polynomials  and  writing 


(5.8) 


(5.9) 


F(<r)  s  (l-crJ^P^)  +  ( 1  +  <r ) n(n  + 1) Pn( t|i)  =0  (5.10) 

where 


*  =  - 


!  1  +<r)b 


[4a  -c  (!+*)*] 


(5.11) 


we  shall  now  show  that  the  roots  of  ( 5. 10)  occur  at  real  values  of  «r  . 
We  observe  that  physically  this  result  is  to  be  expected  because  it  means 
that  the  periods  of  free  oscillation  of  the  fluid  in  the  spheroidal  envelope 
are  real.  In  order  to  prove  this  result,  we  note  that  (5. 10)  is  effectively 
a  polynomial  of  degree  n  +  1  in  <r  and  hence  it  is  enough  to  prove 
that  F(<r)  has  n  +  1  zeros  on  the  real  axis  of  <r  .  One  such  zero  is 
clearly  at  o-  =  -1  ,  where  =  0  .  Further,  between  <r  =  -1  and  <r  =  +1, 
where  ^  1  >  Pn  has  n/2  zeros  and  consequently,  from  the  interlacing 
property  of  the  zeros  of  P  and  P'  ,  F(o-)  has  4  n  -  1  zeros  between 

the  first  positive  zero  of  P^  ,  qua  function  of  and  tr  =  1  .  A  similar 
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remark  applies  to  -3  <  <r  <  -1  .  Again  if  n/2  is  even,  F (tr )  >0  if 
0-+  1  is  small  and  positive,  and  is  negative  at  the  first  zero  of  Pn  for 
which  <r>-l.  Hence  F(c)  vanishes  in  this  range  too.  A  similar  argument 
applies  if  n/2  is  odd.  Counting  up  we  see  that  there  are  exactly  n  +  1 
zeros  of  F(<r)  in  -3<<r<l  all  simple,  and  leading  to  finite  oscillations 
of  the  fluid  except  when  they  coincide  with  poles  of  Bn  . 

Now  consider  the  singularities  of  B  There  are  two  simple  poles 
at  s  -o  and  at  s  =  ik  ,  from  (3.8)  and  (4.10).  In  addition  6^  and  6^ 
vanish  within  the  range  -3i  <  s  <  i  for  any  acceptable  value  of  p  .  In  the 
evaluation  of  B^  we  need  the  weighted  integral  of  u^  ;  consequently 
the  associated  singularities  of  Bn  only  occur  at  s  =  -3i  and  s  =  i  ,  and 
will  be  considered  at  the  same  time  as  the  singularities  of  the  governing 
equation  (2.18).  The  poles  of  Bn  at  s  =  0  and  s  =  ik  can  lead  to 
motions  which  are  apparently  large  as  t  -*  »  if  one  of  the  zeros  of  (5.  8) 
occurs  at  the  same  value  of  s.  The  second  possibility  (s  =  ik)  is  unlikely 
if  n  #  2  since  it  implies  that  two  free  periods  of  oscillation  of  the  fluid  are 
equal  but  we  have  not  been  able  to  rule  this  out.  The  first  possibility  (s  =  0) 
is  real,  and  there  seems  little  doubt  that  for  any  n  we  can  choose  a  value  of 
a/b  to  satisfy  (5.  8)  by  s  =  0;  for  example,  if  n  =  4  and 

6a  =  (n/39  +  vl5)b  , 

the  left-hand  side  of  (5.  8)  vanishes  at  s  =  0  .  Using  the  argument  above 
2  2  2  2 

for  s  -*  ik  ,  (k  =  ( a  -  b  )  /  ( a  +  b  ) ,  we  can  say  that  if  a  =  b  it  is  most 
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unlikely  that  (5.8)  is  satisfied  by  s  =  0  for  n*2.  In  the  case  of 
particular  Interest  therefore,  when  a«b,  we  conclude  that  the  possibility 
of  satisfying  (5. 8)  at  a  ■  0  is  not  likely  to  be  serious  except  if  n  =  2  . 

Let  us  consider  the  case  n  =  2  in  some  detail.  Here  A2  has  a 
double  pole  at  s  =  ik  because  the  coefficient  of  A2  vanishes  and  because 
B2  has  a  simple  pole  from  (4.10).  Using  the  formula  for  A  in  (3.5) 


(5.13) 


The  corresponding  value  of  A2  is  now  easily  worked  out  from  (5.  6)  and, 
reverting  to  the  coordinates  (r,  8,  z),  the  corresponding  contribution  to 


3rz(s2  +  l)B2 
a2(s  -  i)  +  b2(s  +  i) 


(5. 14) 


Let  the  coefficient  of  A  in  (5.13)  be  x  .  The  implication  of  ( 5. 4)  is 
that  the  secondary  motion  described  by  (3. 3)  and  (3.  5)  gives  rise  to  a 
complicated  tertiary  motion  of  which  the  second  harmonic  is  of  the  same 
formas  (3.3)  but  with  a  double  pole  at  s=ik.  Although  formally  this 
means  that  the  tertiary  flow  increases  without  limit  as  t  -*  »  it  is  to  be 
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expected  that  the  double  pole  arises  through  an  error  of  order  x  In  the 
position  of  the  simple  pole  of  A  in  (3. 5).  To  see  this  we  note  that 
the  contribution  from  ( 5. 14)  to  the  tertiary  flow  is  ultimately  dominant 
and  it  gives  rise  via  the  boundary  layer  to  a  similar  form  to  (5. 14)  in 
the  quaternary  flow  but  now  with  a  triple  pole  at  s  =  ik  ,  and  so  on. 
Adding  up  all  these  contributions  to  V  which  are  proportional  to  rz  , 
qua  functions  of  r  and  z ,  their  total  is 

t-i  —  im 

m=0  a  (s-i)  4b  (s  +  i) 

_ 2Qa2i  sin  e(s2  +  l)rz _ 

ws[a2(s  - i)  +  b2(s  +  i)  -3(s2  +  l)x] 

Thus  the  effect  of  the  boundary  layer  is  to  shift  the  pole  s  =  ik  to 


...  .2.  2.2 

,ik 

a  +  b  a  +  b 


the  neglected  terms  being  O(v^)  and  most  probably  0(  $~)  from  (9). 
Of  particular  interest  is  the  position  of  the  pole  when 


1  » 


a2  -  b2 
a2  +  b2 


» 


V 

2 

ua 


(5. 15) 


t 


(5.16) 


the  second  condition  being  fundamental  to  the  present  theory.  So  far  as  x 
is  concerned  it  is  sufficient  to  set  ash  when 


5al  r 

Xs-—  J 

-1 


(5.17) 


where 


62  =  £<“>*  (■♦1-2*0* 


and  since  s  **  0 


^(J7-)ieiir/4(l-2n)i,  if  My, 


It  follows  that 


1)*.  if  ,>fc. 


(5. 13) 


X  =  T"  ®2(  ~~ 7  )*  (°*  195  +  1.  976  i)  , 
2wa 


whence  the  pole  is  moved  to 
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8  a 


,  2  .2. 
(a  -b  I 
2  2 
(a  +*>  ) 


(  )*  (3.720  -  0.3651)  +  ... 

2wa 


and,  since  It  Is  in  the  half  plane  <£(a)  <  0  ,  the  associated  contribution 
to  the  flow  dies  out  as  t  -*  *> .  A  similar  result  can  be  expected  for  all 
b/a  . 

The  last  cases  to  be  discussed  are  the  behaviour  of  the  solution  in 

the  neighbourhood  of  s  »  ±i  and  -31,  where  the  differential  equation 

2  2  2 

(2.18)  takes  on  a  singular  form.  Although  vanishes  at  pa  =  (b  +c  p 
when  s«l  it  follows  from  (4.10)  and  (5.7)  that  Br  is  bounded  as 
s  **  i  for  all  n.  Hence,  substituting  into  (5.6)  and  noting  that 
X4=1y4  on  S0  if  s=i,  it  follows  that 

An  =0(8-1)*  , 

near  s  :i.  Further  the  nth  harmonic  in  (5.4)  can  be  written  as 


^[(s  +  i)2  +  4]*(s  +  i)rzAnGn/2  {^[(s  +  i)2  + 4]r2,  ^(s  +  i)2z2}  , 

where  Gn^2(a,  P)  is  a  homogeneous  polynomial  of  degree  n/2  in  a  and 
with  constant  coefficients.  Hence  the  contribution  to  the  tertiary  flow  is 


(5.19) 


l 

i2 

(5.  20) 

(5.  21) 

P 


#374 


-37- 


O(s-i)  as  s  -*  i ,  and  its  inverse  with  respect  to  s  tends  to  zero  as 
t  -*  «e .  a  similar  remark  applies  to  the  neighbourhood  of  s  =  -3i  .  In  the 
neighbourhood  of  s  =  -i,  B  is  0(s  +  i)  from  (4.10),  and  the 
coefficient  of  in  (5. 6)  is  0(1)  since  n  is  even.  Consequently 

3 

(5.21)  is  0(  s  +  i)  near  s  = -i  and  its  contribution  to  the  tertiary  flow 

tends  to  zero  as  t  -*  » . 

Summarizing  we  can  expect  that,  of  all  the  components  of  the 
secondary  flow  which  depend  on  0,  only  the  steady  component  will  remain 
as  t-»  »fthe  boundary  layer  serving  to  damp  out  the  oscillatory  terms. 


-38- 


#374 


6.  THE  TERTIARY  INVISCID  FLOW  (ii) 


In  this  section  we  consider  the  inviscid  flow  engendered  by  (4.17) 
which  is  independent  of  8  .  We  assume  that  the  associated  inviscid  flow 
is  also  independent  of  8  whence,  on  taking  the  Laplace  transform  with 
respect  to  «t ,  the  governing  equations  reduce  to 


su56  +  2u5r  *  0  • 


(6.1) 


those  equations  also  follow  from  (2.12)  on  assuming  that  ^  is  independent 
of  0  and  neglecting  the  initial  and  forcing  terms.  Further,  on  SQ  , 


r 

2 

a 


,  z  * 

+  V*  U5z 

D 


■  l  *  f?  ^  U5X  ,,b>  “>  > 
a  b  1 


from  (4.17). 


In  terms  of 


s2+  4 


s  9z  ’ 


(6.  2) 


(6.3) 


so  that 


2  ...  8* 


±[rZl]  +  K .  o 

8r  lr  8r  J  +  2  a_2 


s  8z 


(6.4) 


Introduce  new  coordinates  fig  ,  X  5  such  that 


+  >  ,*KX.  ’ 


S  +4 


(6.5) 


2  12  2  1 
where  Yg  =  a  +  4  8  c  »  s0  is  given  by  X  g  =  j  bg  and,  on  SQ  , 


**5  =  *4  a  **1  =  * 


The  boundary  conditions  now  reduce  to 


1  ,.2  2  lA  .... 

;<b  +e  >  u5;  • 


(6.6) 


on  S_  .  The  most  general  appropriate  solution  for  V.  in  S.  is 

0  D  U 


V‘5  *  2„°„P„<  VW  - 


(6.7) 


where  C  are  constants  and,  since  n  *  is  an  even  function  of  p, 
n  ’  5X 


n  also  must  be  even. .  Hence 
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2  n  pi  i  5 .  _  Oicosa 
(2n  +  l)  Vn'2v5'  "  '  2o>s 


(6.8) 


The  main  interest  here  centers  on  the  behaviour  of  the  solution  at  large 

times,  contributions  to  which  arise  from  three  sources.  To  begin  with  C 
’  n 

will  develop  a  simple  pole  whenever  P"  vanishes  unless  s  =  0  :  the 

n 

corresponding  motion,  according  to  (6.7),  is  a  feeble  oscillation  and 
one  can  expect  to  show,  using  a  similar  argument  to  that  which  led  to 
(5.19),  that  viscosity  damps  it  out  ultimately.  Further  the  differential 
equation  is  singular  at  s  =±2i  where  ibs/2yg  =±1  •  Since  P^(±l)  **  0  > 
the  leading  terms  from  this  cause  must  become  zero,  as  t  ■*  «> .  Finally, 
one  must  consider  the  case  s  =  0  .  When  s  is  very  small 


C 

n 


t  Zn  + 1)  flcos  a  (  v  Ibf+cVji  _ 

n(n  +  l)bws  Pn(0)  wa  -1  (-2inar 


(6.9) 


e 

so  that,  on  inverting  (6.  7),  V.  contains  a  term  which  is  proportional  to 
t  when  t  is  large.  The  implication  is  that  the  boundary  layer  exerts  a 
decisive  influence  on  the  inviscld  flow  outside  in  this  instance. 

The  reason  for  the  difference  between  the  solution  in  §5  near 
s  =  -i  and  the  solution  in  §6  near  s  =  0  ,  although  both  governing 
equations  are  similar  and  singular  at  these  points  is  two-fold.  First 
(4. 10)  contains  a  factor  (s  +  i)  while  (4.17)  contains  a  factor  s’^ 
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9 

to  that  tha  bahaviour  of  Vg  near  s  ■  0  is  bound  to  ba  more  singular 
than  tha  bahaviour  of  V*  near  a  «  -i .  There  is,  however,  a  second 
deeper  reason.  Near  s  ■  0  the  governing  equation  (6,1)  reduces  to 


“S*r 


,u& 

02 


(6. 10) 


which  means  that  u  *  is  independent  of  z  and  therefore  must  be 

oz 

an  even  function  of  »*•  Consequently  V5  is  an  odd  function  of 
l*  but  the  boundary  condition  ( 6.  2)  requires  it  to  be  an  even 
function.  This  patent  contradiction  means  that  a  serious  breakdown 

in  the  solution  must  occur  as  t  —  * .  From  studies  of  the  steady  state  (e.g. 
Proudman  (8) )  we  know  that  the  boundary  layer  in  such  a  case  exerts  a 
decisive  Influence  on  the  flow  outside,  adjusting  it  until  the  conditon  that 
u^  is  an  odd  function  of  p  is  satisfied,  which  is  only  possible  if 
u5^  ■*  0  as  t  -*  oo.  Hence  the  initial  angular  velocity  u0  =  Qr  cos  a 
is  also  damped  out  by  the  boundary  layer. 

On  the  other  hand  near  s  *  -i ,  while  the  governing  equations 
(2. 17)  and  (2.18)  also  reduce  in  part  to 


8u  .* 
_ 4z 

02 


=  0  , 


(6.11) 


so  that  u  A  is  an  odd  function  of  n  ,  this  condition  is  actually  satisfied 
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by  (4. 10)  so  that  no  corresponding  difficulties  occur. 
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